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Abstract. 

The Ising-like anisotropy parameter S in the Kondo necklace model is analyzed using 
the bond-operator method at zero and finite temperatures for arbitrary d dimensions. 
A decoupling scheme on the double time Green's functions is used to find the dispersion 
relation for the excitations of the system. At zero temperature and in the paramagnetic 
side of the phase diagram, we determine the spin gap exponent vz « 0.5 in three 
dimensions and anisotropy between < 5 < 1, a result consistent with the dynamic 
exponent z = 1 for the Gaussian character of the bond-operator treatment. At low 
but finite temperatures, in the antiferromagnetic phase, the line of Neel transitions is 
calculated for (5 <C 1 and (5 ~ 1. For c? > 2 it is only re-normalized by the anisotropy 
parameter and varies with the distance to the quantum critical point QCP \g\ as, 
T/v oc l^l^^ where the shift exponent i/j = l/{d— 1). Nevertheless, in two dimensions, 
long range magnetic order occurs only at T = for any S. In the paramagnetic phase, 
we find a power law temperature dependence on the specific heat at the quantum 
liquid trajectory J/t = {J/t)c, T 0. It behaves as Cy oc T'' for (5 < 1 and J sa 1, in 
concordance with the scaling theory for z = I. 
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1. Introduction 

Quantum phase transitions (QPT) from an antiferromagnetic AF ordered state to 
a nonmagnetic Fermi liquid (NFL) in heavy fermion (HF) systems have received 
considerable attention from both theoretical [1] and experimental points of view[2]. 
In contrast to classical phase transitions (CPT), driven by temperature, QPT can 
be driven by tuning an independent-temperature control parameter (magnetic field, 
external pressure, or doping). The physics of HF is mainly due to the competition 
of two main effects: the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction between 
the magnetic ions which favors long-range magnetic order and the Kondo effect which 
tends to screen the local moments and produce a nonmagnetic ground state. These 
effects are contained in the Kondo lattice model (KLM) Hamiltonian in which, only spin 
degrees of freedom are considered. Here we investigated a simplified version, the so-called 
Kondo necklace model [3] (KNM) which for all purposes can be considered yield results 
similar to the original model. While the ground state properties of this model has been 
investigated rather extensively, by a variety of methods H [5l[6l[71[8l[9l[l0l[ni[T2l[T3l[Tl], 
thermodynamic and finite temperature critical properties, close to a magnetic instability, 
remain an open issue. That was warned for us, and it was our first motivation for 
studying the quantum critical properties of this model, as a function of the distance 
to the quantum critical point l^^l at zero and low temperatures [H [15]. We extend now 
this treatment to finite inter-site anisotropy 5 such that, < 5 < 1 since the 5 = 1 
case is appropriate to describe compounds where the ordered magnetic phase has a 
strong Ising component. However, the main reason for considering anisotropy 5 in 
the KNM is to try to describe its effects in the neighborhood of a magnetic quantum 
critical point (QCP) in HF systems, rather than a symmetry problem[12]. This is a goal 
in HF systems, and already several theories were formulated to explain their unusual 
properties [lEl [TTl [18]. Besides, in an early work we were succeeded in finding that the 
Neel line exists since turning on a geometric anisotropy [19], stressing that anisotropy 
is an inherent ingredient in real HF systems. Henceforth, the model will be called 
anisotropic Kondo necklace model (AKNM). This model was already investigated using 
the real space renormalization group machinary|20] but just in one dimension and zero 
temperature. We use the bond-operator approach introduced by Sachdev and Bhatt [2T] 
which was employed previously to both, KLM[22J and KNM[Tl] models but always at 
(T, 5) = (0,0). We find that this method yields a shift exponent that characterizes 
the shape of the critical line in the neighborhood of the QCP, as well as, the power 
law temperature dependence on the specific heat along the so-called quantum critical 
trajectory J/t = {J/t)c, T —>■ 0. We consider the following AKNM: 

H = tY.ir^r^ + i^-S)r!rJ) + Jj:S,.T,, (1) 

<i,j> i 

where and Sj are independent sets of spin-1/2 Pauli operators, representing the 
conduction electron spin and localized spin operators, respectively. The sum {i,j) 
denotes summation over the nearest-neighbor sites. The first term mimics electron 
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propagation which strength t and the second term is the magnetic interaction between 
conduction electrons and locahzed spins Sj via the Kondo exchange couphng J ( J > 0). 
The Ising-hke anisotropy parameter 5 varies from the full anisotropic case 5 = lio the 
well established case 5 = 0. Considering the bond-operator representation for two spins 
S = 1/2, Ti{Si)° = =F|(s|ti,a + tl^^Si ± i(-ap-y't\,i3^i,'y) ('^ = 3^, 2/, -2^) |21j , the Hamiltoniau 
above, at half-filling, i.e., with one conduction electron per site, can be simplified and 
the resulting effective Hamiltonian Hmf with only quadratic operators is sufficient to 
describe exactly the quantum phase transition from the disordered Kondo spin liquid 
to the AF phase, as discussed below. Then, we have a mean-field Hamiltonian: 



k 

+ E [KAjKy + K {tlJ-Ky + ^k,,t-k,,)] , (2) 

k 

where = luo + 2Ak, A'^ = luo + 2A'^, = \ts'X{k), A'^ = \ts^X{k){l - 6) 
and A(k) =X]f=iCosfcs. s is the singlet order parameter consistent with the strong 
coupling limit J/t —>■ oo, where the model becomes trivial, since each S spin captures a 
conduction electron spin to form a singlet, and where the ground state corresponds to 
a direct product of those singlets. The chemical potential /x was introduced to impose 
the constraint condition of single occupancy, is the number of lattice sites and Z is 
the total number of the nearest neighbors on the hyper-cubic lattice. The wavevectors 
k are taken in the first Brillouin zone and the lattice spacing was assumed to be unity. 
This mean-field Hamiltonian can be solved using the Green's functions to obtain the 
thermal averages of the singlet and triplet correlation functions. These are given by. 



{iu^ - uj'^){iu + Ak) 



1 



where ^ = (cu^ — — uj'^). The poles of the Green's functions determine the 

excitation energies of the system as ujq = (-f + /^), which is the dispersionless spectrum 



of the longitudinal spin triplet states, Uk = ±v A| — (2Afc)2 that correspond to the 



excitation spectrum of the x-transverse spin triplet states and uj'f^ = '^\J A'^ — (2A'^)2 
that correspond to the y-transverse one. 
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2. Paramagnetic State 



From these modes above and their bosonic character an expression for the paramagnetic 
internal energy at finite temperatures can be easily obtainedjll [151 [T9| - 

U = Eq + Y^ {uQn{uQ) + t^k^(t^k) + tUk'^(t^k)) (4) 

k 

where Eq = N | Js^ + yus^ — jJ^ +Z]k(^k+t^k~^k~A|^)/2 is the paramagnetic ground 

state energy, n{uj) = \ {coih.^ — the Bose factor, (3 = l/ksT, ks the Boltzman's 
constant and T the temperature. After some straightforward algebra[Tl [19] using Eq. 
(jlj), the paramagnetic free energy renders 

F = 50 - E ln[l + nM] - i E + ^K)] - J ln[l + n^)]- (5) 

For obtaining and fi we minimize the free energy by the saddle-point equations 

2(2 -t) - ^ E coth ^ + coth ^ j + /(.o), 

^ - ^ E f -A(k) coth ^ + ^A(k)(l - 5) coth M") , 

(6) 

where /(u;o) = f (coth ^ - l). 
^. i. Numerical results at T = 

We first study the case T = e.i., without thermal excitations. At zero temperature 
the self-consistent equations given by Eqs. (jS]) can be simplified as. 



4(2 - s') = h{y) + h{y) + h{y) + h{y) 
AJy /3 /i 



4 J 



h{y)-h{y) + h{y)-h{y): (7) 



with 



7-/\ I d k T / \ ^ ^ ^ 



TT" JO 



'l + yX{k) ^''Jo ^i + ^(i-5)A(k) 

hiy) = r d'^k^l + ym, h{y) = rd''k^l + y{l-6)X{k), (8) 
it"- Jo ^ vr" JO * 

where we have introduced a dimensionless parameter y = ts"^ /ujq. An equation about y 

can then be obtained: 
2t 

y = jii-[h{y) + h{y)]/^). (9) 

We will now obtain the numerical solutions to the zero temperature self-consistent 
equations ([7]) using Eq. ([9]). In this case (paramagnetic phase), we have the 2;-polarized 
branch of excitations has a dispersionless value ujz{k) = ujq and the other two branches 
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show a dispersion which has a minimum at the AF reciprocal vector Q = (vr, tt, tt) in 
three dimensions {3d). The minimum value of the excitations defines 

= iJo^l- yd, Ay = uo^l - yd{l - 6). (10) 

The spin gap energy and define the energy scale for the Kondo singlet phase, for 
< 6 < 1 and 6 < respectively. For 6 = 0, A^ and A^ are identical and we obtain the 
original spin gap in the KNM[Tl]. Although we are interested in the case < 5 < 1, we 
consider 6 < due to theoretical reasons. This case only will be consider at T = and 
will not be sketched in this report. The analysis of the spin gap is important because 
the vanishing of gap and the appearance of soft modes define the transition from the 
disordered Kondo spin liquid to the AF phase at the QCP (J/t = {J/t)c,T = 0). At 
this point, it is suitable to clarify that in figures (1), (2) and (3), we sketched the spin 
gap energy like A/ J versus t/J by following the 6 = case[TT], despite we consider 
throughout this paper the control parameter as J/t. That will not yield any physical 
difference since it only will change the onset of the curves from the left to the right. 
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Figure 1. (Color online) The spin gap A/J vs the control parameter t/J is sketched 
for different values of 6 in one dimension and T = 0. ft shows that spin gap is always 
nonzero for < 5 < 1. 

In the one dimensional (Id) case, the energy gap falls linearly for small values of 
t/J and deviates considerably from the linear behavior as t/J gets larger, as it is relates 
in Fig (1). Thereby, it is always nonzero for any 6, supporting its disordered phase, own 
oild Kondo lattices [TTl [23]. 
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t/J 

Figure 2. (Color online) Sketch at zero temperature of the spin gap vs the strength t/J 
in two dimensions. The inset shows the log-log plot of the spin gap versus \{t/ J)c — t/ J\ 
for < 5 < 1. It shows that A/ J vanishes close to {t/ J)c with a exponent vz \. 

The anisotropy dependence on the spin gap in two dimensions {2d) is sketched 
in Fig. (2). For < 5 < 1 the effect of anisotropy is still weak and it changes the 
QCP slightly, until 5 = 0, where both soft modes, and A^, contribute and the QCP 
undergoes a slight jump. Then, the qualitative behavior is the same for this range and 
the gap exponent is approximately z/z ~ 1. It is plotted in the inset of Fig. (2). On the 
other hand, for 5 < 0, the QCP is reduced and the Kondo spin liquid phase is limited 
to a narrower region. It is not shown in Fig. (2). 

In three dimensions, the effect of anisotropy on the spin gap is similar as in the 
2d case. The spin gap follows a exponent vz ~ 0.5 for < 5 < 1 and changes its 
universality for 5 < 0, where the spin gap vanishes close to the QCP more faster. As in 
the 2d case, there is a jump for 5 = which is the result of all soft modes coincide. 

We conclude that, for all anisotropy between < 5 < 1, there exists a critical 
value (t/J)c, where the spin gap vanishes as A/ J oc \(t/J)c — t/J\'^^, and a QPT to 
the ordered magnetic phase occurs in 2d and 3d whereas no transition happens in Id. 
This is similar to the results in Ref. [11] for 5 = and it gives us a kind of universality 
similar as in the isotropic Kondo lattices[Il III [22]. From relation between the spin gap 
and the distance to the QCP, sketched in the onset of Fig. (3), it is shown that when 
t/J increases from its strong coupling limit, the triplet spin gap at the wave vector 
Q = (7r,7r,7r) decreases and vanishes at t/J = (t/J)^. Since A/ J oc \(t/J)c — t/J\^-^, 
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t/J 

Figure 3. (Color online) Anisotropy dependence on the spin gap vs the strength t/J 
in three dimensions at T = and < (5 < 1. The scaling of gap close to the QCP is 
shown in the inset of the figure. It shows the log-log plot of A/ J vs \{t/ J)c — t/ J\ for 
5 = 0,0.1,0.4,1, and scales close to the QCP like A/J ~ |(i/J)c - i/Jj"^ with spin 
gap exponent vz k, 0.5. 

close to the QPT, we can immediately identify the spin gap exponent vz ~ 0.5 at the 
QCP of the Kondo lattice, confirming our early theoretical results[I]. Finally, for 5 < 
exists also a QPT in ci = 2, 3 but no phase transition appears in \d. 

2.2. Analytical results at the quantum critical trajectory 

Since quantum phase transitions are generally associated with soft modes at the 
QCP, where the gap for excitation vanishes, then physical quantities have power law 
temperature dependencies determined by the quantum critical exponents [2¥j: one of 
them is the specific heat Cy, that we will calculate here. This strategy has been 
intensively explored in the study of heavy fermion materials, in the so-called quantum 
critical trajectory J/t = {J/t)c, T — 0, fixing the pressure (in our case the control 
parameter J/t) at its critical value for the disappearance of magnetic order [25]. Then, 
we calculate analytically, the anisotropy dependence on the specific heat at J/t = {J/t)c, 
T ^ for both cases, 6 <^ 1 and 5 ~ 1. All the calculations will be done considering two 
essential approximations: (i) The system is at the quantum critical point J/t = {J/t)c, 
and temperatures T ^ 0. (ii) The temperatures region where will be found the specific 
heat will be lower than the Kondo temperature (Tk). We will begin writing k = Q + q 
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and expanding for small q: \{q) = —d + + O(g^), this yields the spectrum of 
transverse spin triplet excitations as, 

^ ^o\/l + l/A(g) = v^A' + Dq\ 

^ + y\{q){l -5) = + D{1 - 5)q^ + ul5, (11) 

where A = is the spin gap energy given by Eq. (fTOj) since 0<(5< 1, D = ujQ/2d 
the spin- wave stiffness at T = 0, and is the z-polarized dispersionless branch of 
excitations. Considering A = 0, at the QCP[2l] in the spectrum excitations Eq. ffTTj) . 
and using Cy = -Td'^F/dT^ in Eq. ^ we get 

= iSv [ ''™"'<"' + '^?>'^'"''" ^ + """" <i2' 

where Sd is the solid angle. Equation f|T2|) yields the expression for the anisotropic 
dependence on specific heat at the quantum critical trajectory, as an contribution of 
bosons tx and ty. 

Case < (5 ^ 1 — Having shown the relationship between the specific heat Cy and 
5 we now discuss the case 5 <C 1. Making change of variables in Eq. (fT2l) we obtain, 

Cy{5 « 1) = ^^t^ [Ti(rf) + ^(T,(d) - 2T,{d)% (13) 

where Ti(rf) = rfxx'^+i sinh-2(x/2), T2(rf) = /o°° rfa;x'='+2 coth(x/2) sinh-=^(x/2) and 
X = (3ujQq/\fZ. In two dimensions we found Ti(2) = 24^(3) and T2(2) = 96^(3), 
where C, is the Riemann zeta-f unction. In three dimensions Ti(3) = 167r^/15 and 
T2(3) = 1671^^/3. For 5 = 0, the spectrum excitations given by Eq. (ITTl) coincide and we 
recover the exact value as obtained in an previous work for the isotropic KNM[T]. 

Case 5 ~ 1 — Here, it is sufficient to consider, ^ = 1 — 5 <^ 1, where ^ is a 
dimensionless parameter that controls the Ising-like anisotropy in this case. Thereby, 
working in analogy with the preceding case, we obtain 

C.(i«l) = ^l||^(^)\.M(2-i). (14) 

where we have already replaced the ^ expression. The results above show that the 
specific heat of the AKNM for 5 <^ 1 and 5 1 is only re-normalized by the anisotropy, 
concluding that Cy oc T'^ at the quantum critical trajectory for 5 ^ 1 and 5^1. 
Notice that this is consistent with the general scaling result Cy oc T'^/^ with the 
dynamic exponent taking the value z = l[2l]. Since 2; = 1, in three dimensions 
deff = d + z = dc = 4: where dc is the upper critical dimension for the magnetic 
transition [24] . Consequently, the present approach yields the correct description of the 
quantum critical point of the Kondo lattices for d > 3. 

3. Antiferromagnetic Phase 

The mean-field approach can be extended to the AF phase assuming the condensation 
in the x component of the spin triplet like: tk,x = VN't5]i,Q + where t is its mean 
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value in the ground state and ?7k,a; represents the fluctuations. Making the same steps 
as before, the internal energy renders 

U' = 4 + XI (^o'^(^o) + c^k^(t^k) + t^W(^k)) 5 (15) 



where = 



I + - + + ^ - \tZ^^) t^] + Ek(wk + cui, - Ak - K)/2 is 
the AF ground state. The free energy is now 

F' = 4 - 1 X: + ^(^k)] - ^ E ln[l + n{u'^)] - J ln[l + n{uoMlQ) 

Minimizing the free energy Eq. (fT6il . using (OF' /dfi, dF' /ds, dF' /dt) = (0, 0, 0), we can 
easily get the following saddle-point equations, 

-2^1, J fM 




1 + ^(1 + 2nM) + ,1 + ^AO^d + 2.(.0) 




^ = Izts^ - J/4, (17) 
with the excitations spectrum of the x-transverse and y-transverse spin triplet states 



given now by, = iZtsyi + 2A(k)/Z and u'^ = iZt^y 1 + 2A(k)(l - 5)/Z, 
respectively. Generally the equations for s and t in Eq. ( |T7I) should be solved and 
for 6 = the results of Ref. ([1]) are recovered. Here, in the magnetic ordered state, the 
condensation of triplets (singlets) follows from the RKKY interaction (Kondo effect). At 
finite temperatures the condensation of singlets occurs at a temperature scale which, to a 
first approximation, tracks the exchange J while the energy scale below which the triplet 
excitations condense is given by the critical Neel temperature (T^r) which is calculated in 
the next section. Thus, the fact that at the mean-field level, both s and t do not vanish 
may be interpreted as the coexistence of Kondo screening and antiferromagnetism in 
the ordered phase [H [TT|, [22] for all values of the ratio J/t < {J/t)^. 

4. Critical line in the AKNM 

Following the discussion above, the critical line giving the finite temperature instability 
of the AF phase for J/t < {J/t)c is obtained making t = 0. Hence, from Eq. (ITTl) we 
can obtain the boundary of the AF state as, 

where g = \{J/t)c — {J/t)\ measures the distance to the QCP. The latter is 
given by, (J/t), = Z[l - 4^Ek(y|^i^ + 7T^3I(fep)]' separates an 
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antiferromagnetic long range ordered phase from a gapped spin hquid phase. Performing 
the same analysis as in sub-section fl2.2p . expanding the spectrum excitations close to 
Q = (7r,7r,7r), Eq. ( fT8|) becomes 

where we have considered that for temperatures ksT -C ojq, f{uJo) goes to zero faster 
than the first term of Eq. ([T5]) . This equation above allow us to obtain the critical line 
in the AKNM as a function of the anisotropy parameter 6. 



4.1. Case < (5 < 1 

We now demonstrate analytically the appearance of a finite Neel line temperature when 
a small degree of anisotropy 6 in y-component spin is turned on. Then, solving Eq. (|T9|) 
for < (5 <^ 1, we get 



\g\ S,Z^/' (ksTV-'l^ ... 5 



$i(ci) + -($2(ci) + 2<l>i(d)) 



(20) 



where $i(rf) = f^ dxx'^-^ (coihl - l) and ^2{d) = dxx'^+^ smh-^{x/2). We 



notice that the integrals $i(c?) and $2(0^) diverge for d < 3 showing that there is no 
critical line in two dimensions at finite temperatures [H |15] for any anisotropy 5 <^ 1, 
in agreement with the Mermin- Wagner theorem|26]. Nevertheless, for d > 3, the 
integrals are finite and the equation for the critical line shows, (Tisf)s<^i oc \g\'^, with 
= l/{d — 1). If we write the equation for the critical line, f{g,T) = 0, in the form, 
{J/t)ciT) — (J/t)ciO) + VqT^/^ = 0, with vq related to the spin-wave interaction, we 
identify the shift exponent, ip = z/{d + z — 2)[27J, that comparing with gives us the 
dynamic exponent z = 1, a. Gaussian result, since the critical line only exists for d > 2. 
The temperature dependence of the function / arising from the spin-wave interactions 
can modify the temperature dependence of the physical properties, as the specific heat, 
at J/t = (J/t)c- However, in the limit T — we can easily see that the purely Gaussian 
results for the specific heat calculated in section (12.21) is dominant, in concordance with 
the mean-field treatment used here. For 5 = 0, we obtain the well established result for 
the critical line in the KNM[T], this due to the fact that the spectrum energy of the two 
excitations coincide. For 6^1, following the same steps as before, it is straightforward 
to show that the dominant is also T/v(5 ~ 1) oc Ig]"^ for d> 3, and no critical line exists 
for d = 2. 

In summary, we have obtained analytically the expression for the Neel line, below 
which the triplet excitations condense, close to the QCP for < 5 <^ 1 and 6^1. We 
have shown that this line does not exist for d = 2 for any value of the anisotropy, as 
we expected, whereas for d > 3, the power dependence on \g\ of the critical line in the 
presence of the anisotropy is the same of the KNM original. Therefore, the criticality 
close the QCP is governed by the same critical exponents of the isotropic 6 = case 
that we have calculated before [1]. 
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5. Conclusions 

In conclusion, we have examined the phase diagram of the Kondo necklace model in the 
presence of an Ising-like anisotropy at zero and low temperatures by means of analytical 
and numerical techniques. At zero temperature we have derived and solved the self- 
consistent equations on the Kondo spin liquid phase for any value of 5. This allows 
us to calculate the anisotropy dependence on the spin gap for d = 1,2,3. In the Id 
case, there is no indication at all suggesting a critical value for t/ J where the gap would 
vanish for any value of anisotropy 5. For d = 2, 3 we found that the anisotropy in the 
range < 5 < 1 dislocates lightly the QCP until 5 = 0, where the spectrum excitations 
coincide. In this range < 5 < 1 the spin gap exponent is approximately the same, while 
for 5 < a like-jump occurs and it belongs to other universality class. In particular, in 
three dimensions, the triplet spin gap for anisotropy < 5 < 1, close to the wave vector 
Q — (tt, tt, tt), decreases and vanishes at t/ J — {t/ J)c with spin gap exponent vz ~ 0.5, 
consistent with the dynamic exponent z — 1 and correlation length u — 1/2, a result 
in agreement with the mean-field or Gaussian character of the approximations we have 
used to deal with the bond-operator Hamiltonian. On the other hand, at low but finite 
temperatures, wc find that in general the dependence on \g\ of the critical line for the 
AKNM, in the presence of the anisotropy, is the same of the KNM original. This implies 
that the critical exponents controlling the transition close to the QCP, for nonzero 5, 
are the same as those of the isotropic case. We have also obtained the thermodynamic 
behavior of the specific heat along the quantum critical trajectory J/t = {J/t)c, T — > 0. 
It has a power law temperature dependence as Cy oc T'^, a result consistent with the 
scaling theory with the dynamic exponent z — 1. Therefore, the most essential features 
of the Kondo lattices, i.e., the competition between a long-range-ordered state and a 
disordered state is clearly retained in the model for < 5 < 1. The qualitative features 
regarding the stability of the AF phase are well displayed in the model and it allows a 
simple physical interpretation of the phase diagram in anisotropic Kondo lattices. It will 
be left to a further work to compare our theoretical results obtained for the AKNM with 
experimental data in order to clarify to what extent the estimates of 6 from measured 
quantities depend on the theoretical tools used. 
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